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Abstract 



We show that skyrmions arising from compact five dimensional models have stable sizes. We 
numerically obtain the skyrmion configurations and calculate their size and energy. Although 
their size strongly depends on the magnitude of localized kinetic-terms, their energy is quite 
model-independent ranging between 50 — 65 times F%/m p , where F n is the Goldstone decay 
constant and m p the lowest Kaluza-Klein mass. These skyrmion configurations interpolate 
between small 4D YM instantons and 4D skyrmions made of Goldstones and a massive 
vector boson. Contrary to the original 4D skyrmion and previous 5D extensions, these 
configurations have sizes larger than the inverse of the cut-off scale and therefore they are 
trustable within our effective 5D approach. Such solitonic particles can have interesting 
phenomenological consequences as they carry a conserved topological charge analogous to 
baryon number. 



1 Introduction 



Non-linear cr-models can have topological stable configurations known as skyrmions [1] . For 
Goldstone fields parametrizing the coset G/H, skyrmions can exist always that tts(G/H) ^ 0. 
Nevertheless, the existence of non-singular configurations cannot be determined within the 
non-linear cr-model, since their size strongly depends on the UV-completion of the model at 
scales around 4:7iF w , where F n is the Goldstone decay constant. 

Five dimensional gauge theories provide UV-completions of the non-linear cr-models. 
Any compact 5D gauge theory, in which the bulk gauge symmetry G is broken down to H 
at one boundary and to nothing at the other, is described at low-energies by a 4D non- 
linear G/H cr-model with F w ~ ^M 5 /L, where M 5 is the inverse squared of the 5D gauge 
coupling and L is the compactification scale in conformal coordinates. The cut-off scale of 
the 5D gauge theories is estimated to be A 5 ~ 2Air 3 M 5 . Therefore, provided that M 5 > 1/L, 
these theories are valid up to energies much larger than the scale 4nF n . Being this the case, 
compact five dimensional gauge theories allows us to address the question of the stability of 
the skyrmion configurations and calculate their properties. This is the subject of this paper. 

We will calculate numerically the skyrmion configuration arising from compact extra 
dimensional models, concentrating in flat and AdS spaces. We will show that these configu- 
rations have stable sizes, calculable within the regime of validity of our effective 5D theory. 
The size, however, depends not only on M 5 and L, but also on the coefficient of the low- 
est higher-dimensional operator of the 5D theory. The energy of the skyrmion will also be 
calculated and shown to be predicted in a narrow range. 

Attemps to address the stability of the skyrmion in extra dimensional models can be 
found recently in the literature [2-6]. Nevertheless, we consider that none of them address 
fully consistently the stability of the skyrmion configuration (to determine the size within 
the higher-dimensional effective theory approach). Our analysis is also the first to exactly 
calculate the skyrmion configuration in 5D models that requires the use of numerical meth- 
ods. 

The interest in the existence of these topological objects is not only theoretical but 
also phenomenological. In the recent years there has been a lot of activity in using extra 
dimensions in order to break the electroweak symmetry. In these examples the existence 
of topologically stable particles can lead to important cosmological implications. Another 

1 Strictly speaking, 5D theories provide only UV-extensions of the non- linear a- models since they are not 
well-defined theories at arbitrarily high energies. 
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interest in skyrmion configuration arises in the context of warped 5D spaces. The AdS/CFT 
correspondence tells us that skyrmions in 5D are the dual of the "baryons" of 4D strongly 
coupled theories. Therefore studying these 5D solitons can be useful to learn about the 
properties of the baryons. 

2 Skyrmions in compact and warped 5D spaces 

We will be looking for skyrmion configurations arising from the global symmetry breaking 
pattern SU (2)^ x SU(2) R — > SU (2)y. These configurations will still be valid for any breaking 
G — > H always that G and H contains respectively SU(2)l x SU(2) r and SU(2)y as a 
subgroup. This includes, for example, the case SU(N)l x SU(N) r — > SU(N)y. 

The class of 5D theories that we will be considering corresponds to the following one. 
They are SU (2) L x SU (2) R gauge theories with metric ds 2 = a(z) 2 (dx^dx^ — dz 2 ), where x 11 
represent the usual 4 coordinates (4D indeces are raised and lowered with the flat Minkowski 
metric) and z, which runs in the interval [z uv , z m ], denotes the extra dimension. We will take 
a(z) > a(z m ) = 1. We will denote respectively L M and Rm, where M = (/i, 5), the SU(2) L 
and SU(2) R gauge connections parametrized by L M = L a M a a /2 and R M = R a M a a /2 where 
a a are the Pauli matrices. The chiral symmetry breaking is imposed on the boundary at 
z = z m (IR-boundary) by requiring the following boundary conditions: 

(Lfj, - Rfj) \ z=Zm = , (L^ 5 + R^) \ z=Zm = , (1) 

where the 5D field strenght is defined as Lmn = 9m Ln — 9nLm — i [Lm > Ln] and analogously 
for R. At the other boundary, the UV-boundary, we impose Dirichlet conditions to all the 
fields: 

L[i |z=z uv — R/j, \ z =z vv = . (2) 

The Kaluza-Klein (KK) and holographic description of the 5D theory described above has 
been extensively studied in the literature [7,8]. From a 4D point of view, the theory resembles 
to large- iV c QCD where the Goldstone bosons or pions are associated to the fifth gauge field 
component, and the p and a\ resonances correspond to the gauge KK-states. At low-energies 
all these 5D theories are described by an (SU(2) L x SU{2) R ) / SU{2) V non-linear cr-model. 

We are looking for static finite energy solutions of the classical equations of motion (EOM) 
of the above 5D theory. 4D Lorentz invariance allows us to take the anstatz L = Rq = 0, 
Ljx = L/i(x, z) and R^ = it^x, z), where and /t label, respectively, the temporal and spacial 
coordinates while x denotes ordinary 3-space. Starting from the usual 5D YM action the 
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energy of this field configuration reads 

/f Zm M 
d 3 x J dz a(z) -j- Tr [L^ + RpRp*] , (3) 

where the indeces are now raised and lowered by the Euclidean 4D metric. Finding solutions 
to the EOM is the same as minimizing the energy functional Eq. (jHJ), which closely resembles 
the Euclidean YM action in 4D. Our problem is therefore very similar to the one of finding 
SU(2) instantons, even though, as we will see later, there are some important differences 
which will not allow us to find an analitic solution. 

The topological charge of the soliton will be defined by 



Q = 32^ / £X I dZ ^ Tl ~ RiX ° R ^ ' W 



which is the difference between the L and R instanton charges. In order to show that Q is 
a topological integer number, and with the aim of making the relation with the skyrmion 
more precise, it is convenient to go to the axial gauge L 5 = R 5 = 0. The latter can be 
easily reached, starting from a generic gauge field configuration, by means of a Wilson-line 
transformation. In the axial gauge both boundary conditions Eqs. (JT|) and (j2J) cannot be 
simultaneously satisfied. Let us then keep Eq. flTJ but modify the UV-boundary condition 
to 

Li \ z=Zvv = % U(x)diU(xy , Ri \ z=Zvv = , (5) 

where Lj and Ri are the gauge fields in the axial gauge and i runs over the 3 ordinary space 
coordinates. The field U(x) in the equation above precisely corresponds to the Goldstone 
field in the 4D interpretation [9] once a static Ansatz is taken. By using a form notation [10] 
A = —iAf l dx' 1 and remembering that F A F = duj^(A), the 4D integral in Eq. (J4|) can be 
rewritten as an integral of the third Chern-Simons form 003(A) on the 3D boundary of the 
space: 

Q = -L f L 3 (L)-u 3 (R)] . (6) 

The contribution to Q coming from the IR-boundary vanishes as the L and R terms in 
Eq. (jBJ) cancel each other due to Eq. (jTJ). This is crucial for Q to be quantized and it is the 
reason why we have to choose the relative minus sign among the L and R instanton charges 
in the definition of Q. At the x 2 —>■ 00 boundary, the contribution to Q also vanishes since in 
the axial gauge 85 Ai = (in order to have F 5i = 0). We are then left with the UV-boundary 
which we can topologically regard as the 3-sphere S3. Therefore, we find 

Q = j WapiHEW*)] = 24^ / d3xeijkTr [ U9 i uiud J uiud kU j ] eZ. (7) 
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The charge Q is equal to the Cartan-Maurer integral invariant for SU(2) which is an integer. 
Then solutions with nonzero Q, if they exist, cannot trivially correspond to a pure gauge 
configuration and they must have positive energy. Moreover, the particles associated to 
solitons with Q = ±1 will be stable given that they have minimal charge. Eq. (|7j) also tells 
us that the non-trivial configuration £/(x) corresponds to a 4D skyrmion with B being the 
baryon number. In a general gauge, the value of U(x) will be given by 



I7(x) 



pjexp -i I dz'R^z') j-pjexp i f dz! L 5 (x, z') \ 

I L J 2UV J J I L J Z UV J J 



where P indicates path ordering. From a KK perspective, the 5D soliton that we are looking 
for can be considered to be a 4D skyrmion made of Goldstone bosons and the massive tower 
of KK gauge bosons. 

We want to find a numerical solution to the 4D EOM arising from the energy functional 
of Eq.([3]) with suitable boundary conditions enforcing Q = 1. To do that, as we will now 
discuss, the axial gauge is not an appropriate choice. We will then restart with in a gauge- 
independent way and specify later on the new gauge-fixing condition. It is expected on 
general grounds that the solution will be maximally symmetric, i.e. invariant under all 
symmetry transformations which are compatible with the boundary conditions. We impose, 
first of all, invariance under "cylindrical" transformations, i.e. combined SU(2) gauge and 
3D spacial rotations. This corresponds to the following Ansatz [11]: 

_ l + <f>%(r,z) <f>i(r,z) , 2 _ \ A^(r,z) 

j 2 £jak%k ~r ^ y 0j a XjX a j ~r ^ •Ej%a 5 

= difef)^, ( 9 ) 

r 

and similarly for R^. We have reduced our original 4D problem to a much simpler 2 di- 
mensional one, being = {r, z} (where r = Vx^) the 2D coordinates. The Ansatz has 
partially fixed the gauge leaving only a U(1)l x U(1)r subgroup. The allowed trans- 
formations are those which preserve the cylindrical symmetry and have the form g^.R = 
exp[i aL,R(r, z)x a a a / {2r)\. Under this symmetry A^' R transform as gauge fields (A^' R — > 
Ap ,R + dpaL,R) and 4>l,r = <fii ,R + *02 i '' R are complex scalars of charge +1. 

In addition to the global symmetries, our problem has also discrete ones. Those are the 
L <-> R interchange and ordinary parity x — > — x. The charge Q is odd under each of these 
transformations, so we cannot impose our solution to be separately invariant under both, 
but only under their combined action. We can then further reduce our Ansatz by imposing 
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Lp,{x,z) = Rp{-x, z): 

A 1 = Af = -A{, A 2 = A* = -A L 2) 

0! = ^ = -^, 02 = 0f = ^. (10) 

Our solution is now fully specified by 4 real 2D functions A^ and 0. We have a residual Z7"(l) 
invariance corresponding to g\ = 9r = exp[i a(r, z)x a a a / (2r)] under which A^ is the gauge 
field and <fi has charge +1. 

By substituting the Ansatz Eqs. ([9]) and ( TlOl) into the energy Eq. ([3]), one finds 

(11) 



/•oo 

E = 16n dr dz M 5 a(z) 

</0 Jz m 



i|^l 2 + i^ + ^(i-H 2 ) : 



In flat space, a(z) = 1, this corresponds to a 2D Abelian Higgs model (with metric g^ v = 
r 2 5^ v ). For a general warp factor we have, however, a non- metrical theory. Substituting the 
Ansatz in the topological charge Eq. (T4]) we find 



2ty 



oo pz m 



Q — — I dr I dze 

'0 Jz m 



da{-i<p*D i> (f> + h.c.) + F i 



fJ,V 



(12) 



The charge can be written, as it should, as an integral over the ID boundary of the 2D 
space. We will choose the boundary conditions in such a way that the first term of Eq. (TT2l 
vanishes, and therefore Q will coincide with the magnetic flux, i.e. the topological charge of 
the Abelian Higgs model. 

In order to solve numerically the EOM associated with Eq. (fTTj) . they must be recast 
in the form of a 2D system of non-linear Elliptic Partial Differential Equations (EPDE). 
The numerical resolution of the 2D EPDE boundary value problem has indeed been widely 
studied and very simple and powerful packages exist. This is the reason why we cannot work 
in the axial gauge, since the EOM for A\ is not elliptic in this gauge. We will impose the 
2D Lorentz gauge condition d^A* 1 = 0. The equations for A p become J u = dp,{r 2 a F^ u ) = 
r 2 a\3A u + dp,(r 2 a)F' lu where □ is the 2D Laplacian and J is the current of the field <fi. In 
this way we have a system of 4 EPDE and 4 real unknown functions which we can determine 
numerically. Since we are counting the two gauge field components as independent functions, 
one can wonder whether the gauge-fixing condition is satisfied. Notice however that by 
taking the dp derivative of the EOM for A p and observing that the current is conserved, 
dyJ v = 0, on the solutions of the EOM for 0, we find an EPDE for d D A u which reads 
[(r 2 a)D + dp i (r 2 a)d pl ](dpA u ) = 0. This equation has a unique solution once the boundary 
conditions are given. If we impose dpA u = at the boundaries, the gauge will then be 
automatically mantained everywhere in the bulk. 
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We will solve the EOM in the rectangle (z G [z vv , z m ] , r G [0, R]) and we will numerically 
take the limit R — > oo. At the three sides z = z IR , z = z vv and r = R we impose the following 
boundary conditions 

01 = [ 4> = -i e i7r2/L 

, r = B: M. = , (13) 

where L is the conformal length of the extra dimension: 

L= dz = z m - z vv . (14) 




Few comments are in order. The conditions on the IR and UV boundary for 0x,2 and A\ 
come respectively from Eqs. and (|2J), while the condition for A 2 arises from the gauge 
fixing. On the boundary at r = R we have chosen the fields <fi and A 2 with a nontrivial 
profile consistent, however, with E — > at R —>■ oo; the condition for Ai comes again from 
the gauge fixing. Our choice of boundary conditions is such that the charge Eq. (|T2|) receives 
contributions only from the boundary at r = R. The r = boundary of our rectangular 
domain is special, given that the EOM become singular there. The software which we will 
employ (FEMLAB 3.1) permits us to extend the domain up to r = because it never uses 
the EOM on the boundary lines; the equations on that points are provided by the boundary 
conditions. Some special care is required, however, to enforce the program to converge to a 
regular solution. We find that to impose regularity the following conditions are needed: 

" <f>x/r -> A 1 

d 1 A 1 = 

The conditions for 12 and A 2 are extracted from Eq. (Q by requiring the gauge fields to 
be well defined 5D vectors, while the condition for A\ is the gauge fixing. To impose these 
conditions and obtain regular solutions we define the rescaled fields \i,2'- 4>i = r Xi an d 
02 = — 1 + rx2- The boundary condition Eq. (jT5|) reads now ^ = A±, X2 = and A 2 = 0, 
plus the gauge-fixing condition d±Ai = 0. The rescaled fields xi,2, together with Ai j2 , are 
actually the ones used in our numerical computations. 

Any gauge-invariant information about the 5D soliton such as the energy and charge 
densities Eqs. ( ITT]) and (1T21) respectively can be directly extracted with our procedure. It 
would be also interesting, however, to know also what the profile of the skyrmion is. To this 
end we have to go to the axial gauge, as we discussed above. From Eq. (jBJ) we have 



U (x) = exp 



—iaiX l /r / dz'A 2 (r,z') 



exp [i f(r)aiX l /r~\ . (16) 
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We immediatly see that our boundary conditions imply /(0) = and /(oo) = — ir, as it 
should be for a skyrmion in which f(r) undergoes a — it variation when r goes from to oo. 

When applying the above described numerical method, however, we do not find any 
solution for either a flat or warped extra dimension. The 5D theory which we are considering 
does not possess any non-singular soliton and the reason is the following [3,4]. The energy 
Eq. ([3]) of a 5D soliton configuration is bounded from below: 

d 3 x J dza(z)^Tr|e A ^F^F^| > 8tt 2 M 5 \Q\ , (17) 

where we take taken, for simplicity, the Ansatz Lppfaz) = Rp,a(— x, z) = i*/«>(x, z) and in 
the last inequality we have used a(z) > a(z m ) = 1. The first inequality can be saturated 
if Fp,c, is self-dual or antiself-dual, while the the lower bound Stt 2 M 5 \Q\ is approached (in 
warped spaces) when these (anti) self-dual configurations are centered at z = z m and tend to 
zero size (since a(z) has its minimum at z — z IR ). This is the case of a 4D YM instanton 
configuration of infinitesimal size, p — > 0, with the identification of the Euclidean time 
with z, and centered at z = z m . This is clear since the limit p — > is equivalent to keeping 
the instanton size fixed and taking to zero all other dimensional parameters, the curvature 
and the extra dimensional size , a(z) — > 1 and z vv — > — oo respectively. In this limit the 
energy Eq. ([3]) corresponds to the 4D Euclidian action that is minimized by an instanton of 
arbitrary size. This means that the 5D soliton that we are looking for is a singular instanton 
configuration. 

An alternative way to see this is by explicitly calculating the energy of the instanton 
configuration in the limit of small size. This is done in the Appendix. For a flat space we 
obtain 

(18) 

while for the AdS space we have 

(19) 

In both cases we see that the minimum of the energy is achieved for p — > 0. 



E(p) ~ 8tt 2 M 5 



1 + 



1 (P 



4 \L 



E(p) ~ 8n 2 M 5 
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2.1 IR-boundary terms 

The problem of stabilizing the radius of the soliton configuration can be overcome if there is 
a repulsive force on the IR-boundary that pushes the instanton to large sizes. The simplest 
possibility is to introduce an IR-boundary kinetic term for the gauge fields: 



5, 



IK 



[ d 4 x ~ Tr [L^IT + R^BT] . (20) 

J 2, z=z m 
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Figure 1: We plot the bulk energy density Eq. ffTTj) of our numerical solution for flat (left 
panel) and AdS (right panel) space. Both plots are obtained for a/(LM 5 ) = 0.1. 



The presence of this term is expected in these theories since it is unavoidably induced at 
the loop level. An NDA estimate gives a ~ 1 / (167T 2 ) . Furthermore, as we will show in 
section [231 Eq. (1201) is the lowest higher-dimensional operator of our effective 5D theory and 
therefore the most important one beyond Eq. ([3]). 



Using Eqs. Q) and ( ITUl) and the fact that 
boundary energy Eq. (1201 as 



A\ = at z 



E 



IR 



lira 



dr 



(^ 2 ) 2 + ^(l 



we can write the 



(21) 



We can see that this term favors large instanton configurations by substituting the instanton 
Eq. ([36]) in Eq. Q2U). We obtain 

£ir(p) = 6tt 2 -, (22) 
P 

that grows for small p. Therefore we expect that the presence of Eq. (|2ip will guarantee the 
existence of non-singular 5D solitons. These configurations cannot be found analytically and 
for this reason we have to rely on numerical analysis. The effect of Eq. (l2~Ti) corresponds to 
a change of the IR-boundary condition for </> 2 . Instead of that in Eq. ([TBI , we have now 



■2<P2 



a 

z=z m M§ 



(23) 
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Figure 2: Left panel: the thick line represents the soliton energy which we obtained numer- 
ically in the case of flat space. This is well approximated by Eq. (1301 (dashed line) and 
Eq. f[3"B"j) (dotted line) for, respectively, small and large values of a/(LM 5 ). Right panel: 
the same plot for AdS; the approximated curves are provided by Eq. f[3"Tj) (dashed line) and 
Eq. (PJ (dotted line). 



2.2 Numerical results 



To find numerically the 5D soliton configuration we use the software package FEMLAB 
3.1 [12]. This is a numerical package for solving EPDE based on the finite element method. 
We have concentrated on two different spaces: 

Flat: a{z) = 1 , (24) 

AdS: a(z) = — with z vv -> . (25) 

z 

Fig. [1] shows an example of the numerical results that we get; for ajiLM^) = 0.1 we plot 
the bulk energy density of the soliton in flat and AdS space. Notice that the size of the AdS 
soliton is sensibly smaller than the flat one and that the energy density vanishes at z — z vv 
in the case of AdS while it does not in flat space. This is consistent with the holographic 
CFT interpretation of the AdS soliton. Being localized at the IR-boundary, the soliton is a 
composite state of the CFT. 

In Fig. [2] we present the total energy of the soliton configuration E to tai — E + Em for 
different values of a/(M^L). This approximately corresponds to the mass of the soliton 
Mn — E tot ai- Although the energy of the soliton is quite sensitive to a/(M&L) for large 
values of this quantity, this is not the case when E tota i is expressed as a function of the PGB 

2 Here we work at the semi-classical level neglecting corrections due to the quantization of the soliton. 
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decay constant F w and the mass of the lowest KK-state m p that are calculated to be 

? 2M 5 7T / ir 2 a V 1/2 r n 

Fl = -±, m ^_^ 1 + _j for flat space, (26) 



2 \ -I/ 2 

Fl = ^p, m p ~ ^ ( 1 + a J for AdS space , (27) 



~ 4M 5 3tt /„ 9tt^ 

TT 

p AL V ' 32M.L 



where F n is normalized such that the coefficient in front of the kinetic term of the Goldstone 
is F%/4. In Fig. [3] we show the ratio M^m p j F 2 for different values of a. We see that this ratio 
only varies a 20% and stays in the range 50 — 65, resulting in a relatively model-independent 
prediction of the soliton mass. We define the radius po of the configuration as the mean 
radius of the topological charge density, i.e. 



dpi-ipDrf + h.c.) + F nD 



(2? 



For the instanton configuration in a flat uncompactified space po coincides with p. The radius 
Po will go to zero as a — > and will diverge for a —>■ oo. In order to show the behaviour of 
Po in these two limits, we plot in Fig. H] the following combinations: po[Ms/(aL)] 1//2 for both 
flat and AdS space and po[M 5 /(aL 4 )] 1 / 5 for flat space only. We will show below how these 
behaviours can be analytically deduced. 

The numerical results presented above have a simple interpretation in the limit in which 
the boundary term Eq. (1201 dominates or not over the bulk Eq. ([3]). This corresponds to 
the limits a 3> M^L and a M^L respectively. We will discuss each in turn. 

a <C M5L: In this limit the energy of the soliton is dominated by the bulk contribution, 



and therefore we expect that the instanton configuration is a good approximation of the true 
solution. In this case we can find the size of the configuration by minimizing the instanton 
energy with respect to p. For flat space, the minimum of Eqs. (TT51) and (l2"2l corresponds to 



3a 



1/5 



that leads to the total energy 

E "^ 8 -" H, '{ 1+ 4{im) i- (30) 

For AdS space, the minimization of Eqs. ffl9l) and (122]) gives 
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LM; 

Figure 3: The ratio M^m p /F 2 is plotted as a function of a/(LM 5 ) for a flat(dashed line) 
and AdS space. The behaviour for large a/(LM 5 ) is consistent with Eq. (1531 while at small 
a/(LM 5 ) we obtain Eq. ([321). 

These analytical results agree in the small a limit with the numerical ones. From Fig. H] we see 
in fact that po scales as (a/(M 5 L)) 1 ^ 5 and (a/(M 5 L)) 1//2 for flat and AdS space respectively 
in the limit of small a. The energy in this limit is approximately given by the instanton 
energy 8tt 2 M^ independently of the metric or the compactification of the extra dimension as 
can be seen from Fig. [2j Finally, the energy can be rewritten as 

F 2 3tt 3 F 2 

M N ~ 2tt 3 ^ (fiat space), M N ~ £ (AdS space) , (32) 

m p 2 m p 

in accordance with Fig. [3j 

a ^> M5L: In this limit the boundary term Eq. ( 1201) is dominat, and one finds that the 
effective low-energy theory below 1/L corresponds to the following one. It is a 4D theory 
with a non-linearly realized SU (2) l x SU(2)r global symmetry in which the SU (2)y subgroup 
is weakly gauged. The gauge coupling is given by g = 1/ yfa and the mass of the gauge boson 
is then m p = gF n where F n is the Goldstone decay constant defined above. This is valid 
independently of the geometry of the extra dimension. One can find this result by ordinary 
KK reduction that below 1/L leads to massless scalars (L 5 — R 5 ), the Goldstones, plus a KK- 
state of (L M + R p ) with mass ~ a/ M 5 /(La) -see Eqs. ( |26l) and (1271) . A better understanding 
of the large a limit, however, can be obtained by following a "holographic" prescription and 
treating the bulk gauge fields and its value at the IR-boundary as distinct variables. The 
bulk plus the UV-boundary sector has a global SU{2) L x SU{2) R symmetry broken down to 
nothing due to the boundary condition Eq. ([2]). Its 4D spectrum contains the heavy states 
of mass of order 1/L and the Goldstone bosons. On the other hand, the IR-boundary (see 
Eq. ([!]) and Eq. (|20|) ) corresponds to a gauging of only the SU{2) V subgroup of the global 
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Figure 4: Left panel: Po[M 5 / \aL)} 1 / 2 (dashed) and p [M 5 /(aL 4 )} 1 / 5 for the case of flat space. 
Right panel: po[M^/(aL)] 1 ^ 2 in AdS. The limiting behaviour of those quantities is consistent 
with Eqs. (j29|), (j3Tj) and (131). 

SU(2)l x SU(2)r with 1/g 2 = a. Therefore the spectrum of the full 5D theory below 1/L 
consists of the Goldstone bosons parametrizing the coset (SU(2)l x SU(2)r)/SU(2)v with 
decay constant F n plus the SU{2)y gauge bosons that acquire a mass gF w < 1/L. Theories 
like these were considered in the past and it was shown that they have stable non-singular 
skyrmions [13]. The mass of the Q = 1 skyrmion was found to be[f 

Fl 



M, 



N 



59.4 



(33) 



In Fig. [3] one can see that our numerical value of M N tends to Eq. (I3"3"j) for large values of a. 
We also see from Fig. H] that the radius of our 5D soliton scales in the large a limit as 

VaL/M 5 ~ l/m p , 



Po 



(34) 



as expected for a 4D skyrmion. Thus, we can conclude that for a 3> M$L the 5D soliton 
corresponds to a 4D skyrmion made of Goldstones and a massive gauge field. 

Summarizing, by ranging a from to oo, our 5D soliton configuration goes from being 
a small 4D instanton configuration to the 4D skyrmion of Ref. [13]. By increasing a we are 
effectively decreasing the contribution to the soliton of the heaviest gauge KK-states that 
reduces to a single KK in the limit of very large a. 



2.3 Calculability in the 5D effective theory 

We have found that once we include the IR term Eq. (120]) into the action a regular static 
solution to the 5D EOM with Q = 1 exists. In the spirit of effective field theories, however, 

3 We are taking the result from Ref. [13] where for a = m 2 p /(g 2 F%) = 1, m p ~ 770 GeV and F w ~ 93 GeV 
they obtain M N ~ 667 GeV. 



12 



the action also contains an infinite number of 5D higher-dimensional operators which are 
suppressed by inverse powers of the cut-off scale A 5 , and it is of crucial importance to 
establish how much our results are sensitive to such operators. 

The contribution to the soliton energy of the higher-dimensional operators is easily es- 
timated by substituting the soliton configuration into the operators. Since these operators 
have at least one more dimension of energy, their coefficients are suppressed by one more 
power of Ag 1 . On dimensional grounds, therefore, the contribution to the energy of the 
new operators carries on more powers of l/(A 5 p ) where the radius po is the typical size of 
the soliton. We must check that l/(Aspo) C 1 in such a way that a sensible perturbative 
expansion, analogous to the E j A5 expantion for a scattering amplitude at energies E < A5, 
can be performed. If this is the case the 5D soliton can be studied with our 5D effective 
theory in the sense that we can compute its properties, to a certain degree of accuracy, by 
only including a finite number of operators into the action. 

By looking at Eqs. (l29l) and (j3l|) . we find that p A 5 ~ (A 5 L) 4/5 and p A 5 ~ (A 5 £) 1/2 
respectively for a flat and AdS space. Therefore, as long as A 5 L <C 1 (as it should for our 
extra-dimensional theory to make sense), we have l/(Aspo) C 1 as needed in order to trust 
our soliton within the 5D effective approach. We can see this in more detail by considering 
a specific 6-dimensional bulk operator such as F MN D R D R F MN that, by the NDA counting, 
is suppressed by M5/A 2 where A 5 ~ 24ir 3 M 5 . The contribution of this operator to the 
soliton energy is of order 5E ~ 87r 2 M5/(A5po) 2 that should be compared with the boundary 
energy Eq. ( J22l . The relative correction 5E/E m {p ) ~ (247r 3 aA5po) _1 plays the role of the 
expansion parameter. Using Eqs. fl29|) and fl3il we get, respectively 

(247r 3 a)" 6/5 (A 5 L)- 4/5 and (24ttVT 3/2 (A 5 Ly 1/2 . (35) 

If a, whose NDA value is l/(167r 2 ), is not unnaturaly small and A 5 L > 1, the contribution to 
the energy of the new operator is safely small. Notice that, contrary to the naive expectation, 
we did not loose too much in our perturbation parameter because of the fact that we needed 
to include the operator Eq. fl20l) to stabilize the soliton. If we would have found a solution 
for a = 0, its size would have been given by po ~ L and we wolud have obtained, instead of 
Eq. (|35l) . a perturbative expansion parameter of order 1/(A^,L). 

We can safely neglect any ci-dimensional bulk and brane-localized operators with, respec- 
tively, d ^ 6 and d ^ 5. The only operator that could be important in our analysis is the 
five-dimensional Chern-Simons (CS) term which, by NDA, is suppressed by M5/A5 ~ a, and 
therefore it can be as significant as Eq. ( 122|) . This term is however absent in our case due to 
its relation with the anomalies of the theory; it vanishes because SU(2) is an anomaly-free 
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group. Even in models where the CS term can be present, like for example in U{2) gauge 
theories, its effect could be as small as that of dimension 6-operators [4], even though the 
situation is not completely clear. We leave the analysis of the CS effects for the future. 

The issue of calculability is particularly compelling in our case, and it is different from 
what happens in the 4D Skyrme model in which all operators give comparable contribu- 
tions. 3 The problem, in the case of the 4D skyrmions, comes from the fact that a higher- 
dimensional operator of the form a (dU) 4 must be added to the Goldstone kinetic term 
F% (dU) 2 in order to obtain a regular solution of the EOM [1,14]. The value of a estimated 
by NDA is F 2 /A 2 , where A ~ 4:7rF n is the cut-off of the chiral lagrangian. We then see that 
F 2 factorizes in front of the action and the cut-off A is the only dimensionful quantity which 
enters into the EOM, setting the radius of the skyrmion to po ~ 1/A. In our 5D case also 
a higher-dimensional operator is needed to stabilize the soliton. Its coefficient is given by 
a ~ M5/A5 and then also the 5D coupling M5 factorizes in front of the action. Nevertheless, 
our soliton configuration not only depends on the 5D cut-off A5, but also on the conformal 
lenght L. This is the crucial difference from the 4D Skyrme model; a combination of L and 
A 5 sets the size of the 5D soliton fixing it to a scale larger than 1/A 5 . 



3 Conclusions 

We have studied skyrmion configurations arising from compact five dimensional models. We 
have shown that the size of these skyrmions is stabilized by the presence of IR-boundary 
kinetic terms. This size is always larger than the inverse of the 5D cut-off scale I/A5, 
and therefore consistent within our 5D effective theory. This is different from previous 5D 
models [3,4] in which the skyrmion size was found to be of order I/A5. We have numerically 
obtained the skyrmion configurations for different values of a (the coefficient of the IR- 
boundary kinetic term), and calculated their size and energy. Although the size of the 
skyrmions depends strongly on a, their energy is quite model- independent. By varying a 
these skyrmion configurations smoothly interpolates between small 4D instantons and 4D 
skyrmions made of Goldstones and a massive gauge boson. 

The existence of stable skyrmion configurations in extra dimensional models rises different 
phenonemological issues. For example, theories of electroweak symmetry breaking arising 
from extra dimensions will have this type of stable configurations, and therefore the analysis 
of their cosmological consequences are of great importance. For 5D Higgsless [7] or composite 

4 When this is the case, not only the effective theory does not give us any quantitative information on the 
soliton, but also its very existence is doubtful. 
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Higgs models [15] in which F n = 246/e GeV (where e < 1) and m p ~ 1.2/e TeV, we find a 
soliton energy Mn ~ (2.5 — 3)/e TeV. These TeV stable particles are possible dark matter 
candidates, therefore the precise determination of their relic abundances is of important 
phenomenological interest. These configurations can also be useful in holographic models 
of 4D strong interactions. The AdS/CFT correspondence tells us that these 5D solitons 
correspond to the "baryons" of the dual theory. The 5D AdS model studied here has already 
been proposed as a holographic model of QCD, giving predictions for the meson spectrum 
and couplings in good agreement with the experimental data [8]. The skyrmion found here 
gives an approximate value for the mass of the proton. We find — 500 — 650 GeV, too 
low compared with the experimental value ~ 1 GeV. We must notice, however, that in our 
analysis we have not included the 5D CS term, responsible in holographic QCD models of 
the Wess-Zumino-Witten (WZW) term. It is known that the presence of the WZW term has 
important impact on the skyrmion mass [16], so we can expect that including the CS in our 
analysis can enhance our prediction of the proton mass. These and other phenomenological 
questions deserve a further analysis that we leave for a future work. 
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Appendix: The instanton energy in compact and warped 
spaces 

In the 2D notation used in this paper the BPST instanton with center at (x = 0, z — z m ) 
and size p corresponds to [17] 

d $ d r ^> 

0! = T^— = 0l(p) , 02 = "I ~r— — = 2 (P), / 

^1 = 1 = ^), A 2 = -^=A 2 (p), 

where 

$ = . (37) 

p 2 + r 2 + (z - z m ) 2 1 ' 

The singular instanton of size p — > solves, up to a gauge transformation, the variational 

problem defined by the energy Eq. (flT|) and the boundary conditions Eqs. (jTBl and (TT5]) . 
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Indeed, it has the same topological charge (Q = 1) of the solution we were looking for and, 
since its 2D energy density is a 5-function at r = 0, z = z m , it has the minimal energy 
E = 87r 2 M5 which saturates the lower bound of Eq. ( fTTl) . 

We can explicitly check that the instanton Eq. (136]) fullfills, for any p, the boundary 
conditions at r = and z = z 1R given in Eqs. ffT3l) and (JT3J) - At the other two boundaries 
the instanton of size p — > has different boundary conditions: 



where j3 = 2 arctan [r/(zi R — z)]. These boundary conditions are, as anticipated, topologi- 
cally equivalent to those in Eq. ( fl3l) and then one can convert one into the other by a gauge 
transformation. Instead of gauge rotating the instanton of p = to make it fulfill Eq. (fi"3l . it 
is more convenient to rephrase our original problem in the new gauge and the new boundary 
conditions Eq. ( |38|) . 

Let us consider instantons of small but non-vanishing size. We would like to compute 
the energy E(p) of such configurations. We already know that E(p) will have an absolute 
minimum at p = 0; this means that a small instanton is subject to an "attractive force" 
which tends to shrink its size to zero. By computing d p E(p) one can measure the strength 
of this force. It is important to remark that there are two different effects which make the 
force arise. The first one, which is essentially local, is due to the curvature of the space and 
therefore it is only present in warped spaces. It pushes the instanton to be localized as much 
as possible at z m where the warp factor has its minimal value and then it is energetically 
favorable. The second effect is non-local and is due to the presence of the boundary at z vv 
which confines the solution into a finite volume. At the more technical level the two effects 
come, respectively, from the fact that the instanton fails to fulfill the bulk EOM when the 
space is warped and the boundary condition Eq. (|38|) at z vv . It cannot then be a stable 
configuration. 

To compute E(p) let us first try to substitute the instanton configuration in the energy 
functional Eq. (FTTl) . We call Eb the "bulk" contribution to the energy which we obtain in 
this way. We will see later that this is not the only contribution to E(p). For a generic warp 




r 



oo 




(38) 
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factor we have □ 

E B {p) = 12 vr 2 M 5 f dy = 12 n 2 M 5 £ p n c n (p/L)^^ , (39) 

where we have expanded a(z) in Taylor series around z = z IR and 

f° y n 

c n (p/L) = dy . . (40) 

J-L/p (i + r) 5/2 

The leading term of Eb in the small p expantion comes from the n = term which reads 

12 n 2 M 5 c (p/L) * 8 tt 2 M 5 (l - |^ + 0(p 6 /L 6 )) , (41) 

where we have also expanded c (p/L) around p = and have kept, for future convenience, 
the first correction of order p 4 . A linear contribution to the energy can only come from the 
n = 1 term and is given by 

-4n 2 M 5 pa^(z m ) (l + 0(p 3 /L 3 )) . (42) 

If is non zero, this term must be negative since the warp factor is a decreasing function 
of z. Therefore the above term gives, as expected, an attractive force which makes the 
instanton shrink. In the case of AdS, Eq. (1231) . we have = —1/L and we obtain the 
result of Eq. (fl9l) . For a space with = the linear contribution vanishes and the leading 
correction to the energy, which comes from the n = 2 term, is of 0(p 2 /L 2 ). Making the space 
less and less warped near z m we are increasing the power of pj L of the instanton energy, and 
therefore making the attractive force weaker and weaker at small p. This process however 
stops at order p 4 /L 4 . The n = 4 term is proportional to C4 ~ In p, and for n > 4 we have 
c n ~ 1/p™ -4 , so all terms n > 4 contribute at the order p A /L 4 to the energy. Therefore we 
expect that the weakest possible force will be obtained from a energy of 0(p A /L A ) that is, 
as we will now show, what actually happens in the case of flat space. 

For flat space only the n = term appears in Eq. ( |39l) and Eb is given by Eq. (T4T1) . We 
then immediatly see that Eb cannot be the total instanton energy. If so, the point p = 
would not be a minimum but a maximum, and we would not find an attractive but a repulsive 
force. The reason why there is an extra contribution to the energy is that the instanton of 
finite size does not respect the boundary condition at z uv . Therefore, it does not belong to 

5 We are considering, everywhere in this Appendix, z uv 7^ because for z uv = the energy of the 
instanton in the AdS slice diverges. Our numerical computations show however that the true solution is 
widely insensitive to the position of the UV-boundary and the results we derive here can be applyied for 
z uv = as well. 
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the set of allowed field configurations on which our variational problem is defined. In order 
to find the instanton energy we need a different formulation of the variational problem in 
which the allowed field configurations do not necessarily respect the boundary condition. 
The latter will arise, like the EOM, from the minimization of a new energy functional which 
is now defined to act on a more general set of configurations, to which the instanton belongs. 
This new functional is given by Eq. (fill) and the UV-boundary term % 

(<j> - m)*D z ct> + (0 - 0(0)) D z <p* A 1 {0)) F l2 . 

(43) 

It is easily understood why the addition of E vv makes the UV boundary conditions arise as 
EOM. The variation of Eq. ( TTTT) and (143j) gives, in addition to the usual bulk terms whose 
cancellation will give rise to the EOM, localized UV-terms of the form (<fi — 0) 8{D z (j)) 
and (Ax — A{j 5(F 12 ). Requiring those terms to vanish enforces the boundary conditions in 
Eq. (EHD. 

We can finally compute the energy E(p) of the finite size instanton. It is given by the 
sum of E B in Eq. fl39l and E vv in Eq. (j4"3l . in which of course we have to substitute the 
instanton configuration. The latter term reads 

E vv (p) = 57r 2 M 5 a(z uv )£ i , (44) 

and gives, for a generic warped metric, a p 4 contribution to the energy. This contribution 
does not affect the result in the case of the AdS slice, since it is subleading, but it is crucial 
in flat space as it corrects the negative sign which we found in Eq. (|4"T1) . Adding Eg and E vv 
we obtain Eq. (fTBl) . 
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